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We study a multiorbital model for the alkaline iron selenides Ki_ a; Fe2- i/ Se2 using a slave-spin method. With 
or without ordered vacancies, we identify a metal-to-Mott-insulator transition at the commensurate filling of six 
3d electrons per iron ion. For Hund's couplings beyond a threshold value, this occurs via an intermediate orbital- 
selective Mott phase, in which the 3d xy orbital is Mott localized while the other 3d orbitals remain itinerant. 
This phase is still stabilized over a range of carrier dopings. Our results lead to an overall phase diagram for the 
alkaline iron selenides, which provides a unified framework to understand the interplay between the strength of 
vacancy order and carrier doping. In this phase diagram, the orbital-selective Mott phase provides a natural link 
between the superconducting Ki_ :r Fe2-! / Se2 and its Mott-insulating parent compound. 



Introduction. An important question of active discussions 
is the strength of the electron correlations in iron-based su- 
perconductors (Hill]. The issue is crucial for understanding 
the properties of both the normal and superconducting states 
of these systems. For the iron pnictides, the antiferromag- 
netically ordered metallic ground state of the parent com- 
pounds J3| may arise either from the Fermi surface nesting of 
a weak coupling theory |4j,|5[], or from the strong correlation 
effects associated with the proximity to a Mott transition and 
the concomitant quasi-local moments J6- 12 1. A number of 
factors argue in favor of the latter picture, including the large 
spectral weight in the fluctuating magnetic spectrum IU3I1 . For 
the 11 iron chalcogenides [14J], both the large ordered mag- 
netic moment and the ordering wave vector are difficult to 
understand within the nesting picture. Also pertinent is the 
iron oxychalcogenide La2C>3Fe2Se2, in which the Mott insu- 
lating behavior has been experimentally identified and theo- 
retically explained in terms of the band narrowing effect as- 
sociated with the expansion of the iron lattice unit c ell 1 1 511 - 
These results suggest that the incipient Mott picture lla 1 1 '711 
is even more pronounced in the iron chalcogenides. 

The recently discovered alkaline iron selenide supercon- 
ductors ll 811 Ai_ x Fe2_ y Se2 (A=K, Rb, Cs, or Tl) shed new 
light on this issue. In these materials the superconducting 
T c is comparable to that of the pnictides Ill8l-| 2 lll. a nd the 
superconductivity is near an insulating phase 11221 12311 . The 
insulator is antiferromagnetic with a very large ordered mo- 
ment 11241 |25J, and is intimately connected to the ordered 
iron vacancies 11221 l24l 12611 . The lack of hole pockets in the 
Fermi surface of the superconducting compounds revealed by 
the ARPES measurements Il27l429ll makes the high T c super- 
conductivity and the large-moment magnetic order hardly ex- 
plainable by the nesting mechanism. Instead, they are more 
naturally understood in terms of the incipient Mott picture. 
For instance, the insulating state is naturally interpreted as 
a Mott insulator (MI), not only because it would have been 
metallic in the absence of interactions but also because the 
interactions are strong as inferred from the large ordered mo- 
ment; the enhanced interaction effects have been attributed to 
the band-narrowing caused by the ordered iron vacancies 1I30I - 
3211 . Various experiments suggest that the superconducting 
state is either free of iron vacancies or vacancy disordered, and 



is intrinsically phase separated from the vacancy ordered insu- 
lating state B3I435I1 . A key open question is how the vacancy 
ordered insulating state connects to superconducting phase. 
Elucidating the linkage between the two phases is an impor- 
tant goal of the present study. 

It is well known that tuning either correlation strength 
or carrier density may induce metal-to-insulator transitions 
(MITs) in correlated electron systems l36ll . In multiorbital 
systems, the physics associated with this transition is richer 
than that of one-band systems since both the bandstructure and 
the correlation strength may be orbital dependent. An extreme 
example is the orbital-selective Mott transition (OSMT), for 
which the Mott transition takes place at different correlation 



strengths for different orbitals [3 
OSMT occurs in (Ca,Sr) 2 Ru0 4 T37 



38h- It is believed that the 
3911 For iron-based super- 



conductors, strong orbital differences have been suggested in 
several systems For the iron pnic ides, it was shown 



theoretically 14311 that the OSMT has nearly as competitive a 
ground-state energy as the other competing phases but is ulti- 
mately not stabilized as a ground state. 

In this letter, we investigate the MIT in the alkaline iron se 
lenides system Ki_.j;Fe2_ySe2 using a slave-spin method 143 
4411 . We show that when the Hund's coupling is sufficiently 



strong, the Mott localization of the system is always via an 
OSMP, in which the 3d xy orbital is Mott localized, while the 
other orbitals are still itinerant. This OSMP generally exists in 
both the iron vacancy ordered and disordered cases, and sur- 
vives a range of carrier doping. It provides a necessary con- 
nection between the vacancy ordered insulating phase and the 
metallic normal state above T c . Our results allow us to make 



contact with recent ARPES measurements in this system 145 1. 

Model and method. We consider a multiorbital Hubbard 
model for the Ki_ a; Fe2_ a Se2 compound. The Hamiltonian 
reads 



H = H + Hi, 



(1) 



Here, Hq contains the tight-binding parameters among the five 
3d-orbitals, 

H = \ *if 4aa d j^ + 5Z( £a ~ M)4«r d ia<7, (2) 

where d\ aa . creates an electron in orbital a = 1, 5 with spin 
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FIG. 1. (Color online) a),c): The regular and 1/5-depleted square lat- 
tices, respectively corresponding to the alkaline iron selenides with 
disordered and the \/E x \/5 ordered iron vacancies. b),d): Corre- 
sponding phase diagrams in the J-U plane at commensurate filling 
N — 6 per Fe for the multiorbital model. The dark and light shaded 
regions respectively refer to the Mott insulator (MI) and the orbital- 
selective Mott phase (OSMP). The red circles and blue diamonds re- 
spectively denote the Mott transition and the crossover between the 
fully itinerant metal phase and OSMP. The black dashed line shows 
the crossover scale U* between the weakly and strongly correlated 
metals. 



a at site i, e a is the on-site energy reflecting the crystal level 
splitting, and \i is the chemical potential. We have taken the 
tight-binding parameters from Ref.|45l and the units are in eV. 
7?i nt contains on-site Coulomb interactions 
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where n JQCr = d\ aa di aa . In this model, U, U', and J re- 
spectively denote the intraorbital repulsion, interorbital re- 
pulsion, and Hund's rule exchange coulping. We will take 
U' = U - 2 J. j46ll 

The MIT of the above model is studied using a U(l) slave- 
spin method 14311 . Here, a slave S=l/2 quantum spin is intro- 
duced to carry the charge degree of freedom, and the metal- 
lic (Mott insulating) state corresponds to the magnetically 
(dis)ordered state of the slave spins with quasiparticle spectral 
weight in each orbital Z a > (Z a = 0). For simplicity, in 
the calculation we drop the spin-flip and pair-hopping terms 
in the interaction H mt . Including these terms leads to simi- 
lar results ll4"oll . We study the MIT on two two-dimensional 
(2D) lattices of iron ions: a regular square lattice sketched in 
Fig. 02 a) and a 1/5-depleted square lattice shown in Fig. [He). 
They respectively stand for the completely disordered and the 
perfect -\/5 x \/5 iron vacancies. 



FIG. 2. (Color online) (a) and (b): Evolution of orbital resolved 
quasiparticle spectral weight Z a (in (a)) and orbital filling factor (per 
iron site per spin) with U for the multiorbital model at N = 6 and 
J/U = 0.2 on the regular lattice. The vertical dashed line indicates 
the position of U* , which specifies the dashed lines in Fig. fTJb)(d). 
(c) and (d): Same as (a) and (b) but on one of the two inequivalent 
sites of a unit cell of the 1/5-depleted lattice. Z a for the xz and yz 
orbitals switch on the other site, as do the orbital filling factors. 



Metal-to-insulator transition. The results at the commen- 
surate filling corresponding to N = 6 3d electrons per Fe 
are summarized in the phase diagrams of Figs. [TJb) and (d). 
In both the vacancy disordered and ordered cases, the system 
experiences a Mott transition at Umt from a metal to a MI 
with increasing U. The insulating phase is a low-spin MI for 
J/U < 0.01, but a S=2 high-spin MI for larger J values. 
Umt first decreases then increases with increasing J/U ratio. 
Such a nonmonotonic behavior is a general feature of systems 
away from one electron per orbital, and is also obtained in the 
five-orbital model for the parent iron pnictides [43]. When the 
Hund's coupling is above a threshold (J/ U > 0.1), the system 
crosses over from a weakly correlated metal to a strongly cor- 
related metal with increasing U. The onset of this crossover 
(at U*) is identified by a rapid drop of Z a and a kink in the 
orbital filling in each orbital, as shown in Figs.|2ja)-(d). In the 
strongly correlated metallic state, Z a is strongly orbital de- 
pendent. Increasing U does not lead to the simultaneous Mott 
localization of all orbitals. The metallic state first crosses over 
to an intermediate OSMP at Uos- The Mott transition then 
takes place between the MI and the OSMP at a larger U. 

Importantly, the phase diagram of the vacancy ordered sys- 
tem is similar to that of its vacancy disordered counterpart. 
Quantitatively, Umt and Uos are respectively smaller in the 
vacancy ordered system, which reflects the ordered-vacancy- 
induced reduction in the kinetic energy and hence enhance- 
ment in the correlation effects 1301. 



Nature of the orbital-selective Mott phase. As shown in 
Figs. 13 a) and (c), in the strongly correlated metallic regime, 
Z in the xy orbital is suppressed the most, and this orbital is 
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very close to half-filling. Further increasing U results in the 
Mott localization of the xy orbital at Uos- The other orbitals 
remain itinerant up to Umt- The system is thus in an OSMP 
for Uos < U < Umt- We now turn to discussing the factors 
that stabilize the OSMP. For simplicity, we limit our discus- 
sion to the vacancy disordered case. The vacancy ordered case 
is qualitatively similar. 

We start from the physics that governs the crossover be- 
tween the weakly and strongly correlated metals. Fig. [2 a) 



plots the effective magnetic moment 5 c ff = y ({S z ) ) as a 
function of J/U. It rapidly increases when the system passes 
through the crossover. In the strongly correlated metallic 
phase, S e s ~ 2, indicating that the 5 = 2 high-spin con- 
figuration, promoted by the Hund's coupling, is dominant in 
this regime. 

The Hund's coupling also suppresses the inter-orbital cor- 
relations C a< p = (n a np) — (n a )(np), as shown in Fig. |3ja). 
Together with the crystal level splitting, this effectively decou- 
ples the xy orbital from the others because in Ki_ x Fe2- y Se2 
the xy orbital is the topmost level and is well separated from 
the others. For the same reason, it is easier to stabilize the xy 
orbital to be at half-filling for an overall filling of six elec- 
trons per Fe, as shown in Fig. |2jb) [and Fig. |2jd) for the 
vacancy-ordered case]. Compared to the degenerate xz/yz or- 
bitals, which are also close to half-filling, the threshold value 
for the Mott localization in the non-degenerate xy orbital is 
smaller. Moreover, in the noninteracting limit the density of 
states (DOS) projected to the xy orbital is narrower than those 
of other orbitals. For instance, as shown in Fig.|3jb), the ratio 
of the width of DOS for the xy orbital to that for the xz/yz 
orbitals is about 0.6. (This is also the case for the vacancy 
ordered model; see the figure in the Supplementary Material.) 
This ratio is smaller than that for the LaOFeAs system, which 
is about 0.7. Hence in Ki_ K Fe2-j,Se2 the xy orbital contains 
less kinetic energy. Taking into account all the three factors, 
we find that in the strongly correlated metallic state, it is much 
easier to drive the xy orbital toward the Mott localization. The 
result is the OSMP. 

The threshold interaction for the OSMP, Uos, shows a 
strong dependence on J especially when J/U is large. This 
seems counterintuitive; one could expect that, if the xy or- 
bital is fully decoupled from the others, Uos should ap- 
proach the critical U of a single-band Hubbard model, and 
hence should not depend on J. To understand the behav- 
ior of Uos, we examine the propagation of a charge exci- 
tation by assuming both U and J are large so that we may 
take the ground-state configuration to be the 5 = 2 high-spin 
state. A charge excitation with one more electron filled in the 
ground state can propagate via hopping to neighboring sites. 
K xy (Kxy) denotes the kinetic energy gain associated with 
the hopping processes (not) involving the xy orbital. Two rep- 
resentative hopping processes that do not disturbe the high- 
spin ground-state configuration are illustrated in Fig. Oc). 
Note that K xy < K^y, not only because the xy orbital 
has a narrower (non-interacting) bandwidth but also because 
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FIG. 3. (Color online) (a): Evolution of the effective moment Ses 
and interorbital correlations C XZtXy and C xz , y z with J/U at N — 6 
and U = 3 eV on the regular lattice. The vertical dashed line indi- 
cates the position of U*(J). (b): Orbital projected density of states 
in xy and xz/yz orbitals at the noninteracting limit ( J=U =0) of the 
same model. The curve for the xy orbital is shifted upward for clar- 
ity, (c): Sketch of the two hopping processes of a charge excitation 
on the high-spin ground-state configuration that avoid penalty from 
the repulsive interactions, (d): The differences in free energy den- 
sity as a function of J/U at U — 4.2 eV, showing the competition 
among the metallic (M), Mott insulating (MI), and orbital-selective 
Mott (OS) solutions. 



the Hund's coupling suppresses the interorbital fluctuations. 
The Mott gap associated with either process is estimated as 
A a = E{N+1)+E(N-1)-2E(N) « U-3J-K a , where 
a = xy,xy. The difference in the kinetic energy gains leads 
to two different Mott gaps. Uos {Umt) can be estimated as 
the U value where the Mott gap A a vanishes. Hence we have 
Uos ~ K xy /{1 - 3J/E7) and U MT ~ K w /(1 - 3 J/17), 
both of which increase with J/U. This general consideration 
is consistent with our calculated phase boundaries for suffi- 
ciently large J/U. Interestingly, in this regime, increasing 
J at a fixed U leads to a derealization from MI to OSMP, 
and then to the metallic phase. This is further confirmed by 
comparing the free energies of the three states, as shown in 
Fig. [3d). Note that in the metallic phase, the inter-orbital cor- 
relations between the xy orbital and others are substantially 
suppressed, but remain nonzero. The above argument on the 
behavior of Uos does not hold for smaller J/U where Uos is 
close to U* . In this regime, since the ground state mixes both 
high- and low-spin configurations, several mechanisms favor- 
ing either increasing or decreasing Uos compete. As a result, 
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FIG. 4. (Color online) (a) and (b): Phase diagrams with U and carrier 
doping concentration 5 for the multiorbital model at J/U = 0.2 on 
the regular (in (a)) and the 1/5-depleted (in (b)) lattices, respectively. 
In either diagram, the large red dot refers to the Mott transition, and 
the blue diamonds shows the orbital-selective Mott transition, (c): 
Sketch of a material-based phase diagram in the plane of carrier dop- 
ing S and ordered vacancies. The vacancy order parameter has been 
scaled to be between and 1 (see text). The vacancy ordered in- 
sulating state is located as the cyan square in this phase diagram. 
The superconducting state is tentatively placed as the orange triangle 
on the 8 axis. The dash-dot line shows a possible route to connect 
the two phases. For realistic parameters, the Mott transition point is 
close to the origin, but could be either above or below it. 



Uos shows complicated, even nonmonotonic, J dependence 
(Fig. Bib)). 

Orbital-selective Mott phase at finite dopings. Unlike the 
MI, which exists only at a commensurate filling, the OSMP 
can be stabilized at incommensurate fillings if the chemical 
potential of the itinerant carriers falls inside the Mott gap of 
the localized orbital. In Figs. |4ja) and (b) we show the U 
vs. carrier doping concentration (6 = N — 6) phase diagrams 
for both the vacancy disordered and ordered cases. In both 
systems, the Mott transition takes place between a MI and an 
OSMP at the commensurate filling 5 = 0, and an OSMT be- 
tween the OSMP and the metal extends to nonzero doping 
concentrations. In the vacancy ordered case, both the Mott 
transition and the OSMT take place at lower U values com- 
pared to the vacancy disordered case, reflecting the enhanced 
correlations due to ordered vacancies. 

Unified phase diagram for alkaline iron selenides. In light 



of the above two phase diagrams, we propose a material-based 
unified phase diagram for both the vacancy ordered and disor- 
dered compounds. In this diagram, sketched in Fig. |4jc), the 
horizontal axis refers to the carrier doping 5, and the vertical 
axis stands for the strength of the vacancy order. In general, 
an Ki_ K Fe2_ y Se2 system contains both vacancy ordered and 
disordered regimes. We may parameterize the strength of the 
vacancy order from to 1, according to the volume fraction of 
the vacancy ordered regime (or, alternatively, to the potential 
strength of a virtual Fe atom, with the vacancy correspond- 
ing to an infinite potential). The two limiting cases along the 
vertical axis in the phase diagram, and 1, are obtained from 
Figs. IHa) an d (b) at a fixed U (which takes the value in real 
materials), respectively. The remaining part of the phase dia- 
gram can then be constructed by interpolating between the re- 
sults in Figs. |4{a) and (b) at the same U. The resulting diagram 
consists of a MI, an OSMP, and a metallic phase. The profile 
is similar to what are shown in Figs. [U a) and (b) because, ef- 
fectively, the ordered vacancies enhance the correlation U/D. 
(Here D is a characteristic bandwidth of the multiorbital sys- 
tem.) The insulating compound with the y/5 x y/E vacancy or- 
der is located at S = and vacancy order 1 (square cyan sym- 
bol in Fig.|4jc)). On the other hand, we tentatively place the 
superconducting phase of the superconducting compounds at 
ambient pressure to be on the S axis in the metallic state close 
to the OSMT (triangular orange symbol in Fig.[4|c)). We see 
that one physical trajectory going from the insulating phase to 
the superconducting one is for the dopants to both introduce 
extra carriers and suppress the vacancy order. During the evo- 
lution of the system, the OSMP is an unavoidable intermediate 
phase connecting the insulatin g an d superconducting states. A 
recent ARPES measurement 14511 has provided evidence that 
the superconducting phase of Ki_j;Fe2-j,Se2 is indeed close 
to the OSMT. 

In summary, we have studied the metal-to-insulator tran- 
sition in the five-orbital Hubbard model for the alkaline iron 
selenides Ki_. E Fe2_j / Se2 with and without ordered vacancies. 
We find that the Mott localization of the system is via an in- 
termediate orbital-selective Mott phase, in which the 3d xy or- 
bital is localized while the other 3d orbitals are itinerant. The 
orbital-selective phase is stabilized by the combined effect of 
the orbital dependence in the bandwidths and the pinning of 
the xy orbital to half-filling due to both the Hund's coupling 
and crystal level splittings. It persists over a range of carrier 
dopings. Finally, we have proposed a unified phase diagram 
for the alkaline iron selenides, with the orbital-selective Mott 
phase providing the link between the insulating and supercon- 
ducting compounds. 
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SUPPLEMENTARY MATERIAL 

Fig.[5]shows the orbital projected density of states (DOS) in xy and xz/yz orbitals at the noninteracting limit (J=U=0) of the 
vacancy ordered model. The DOS has been averaged over the two inequivalent sites of a unit cell. In the vacancy ordered model, 
the DOS in the xy orbital is narrower than those in the other Fe 3d orbitals, and the ratio of the width of DOS for the xy orbital 
to that for the xz/yz orbitals is about 0.6. These are both similar to the vacancy disordered model. 
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FIG. 5. (Color online) Orbital projected density of states (DOS) in xy and xz/yz orbitals at the noninteracting limit (J=U=0) of the vacancy 
ordered model at N = 6. The curve for the xy orbital is shifted upward for clarity. The DOS has been averaged over the two inequivalent sites 
of a unit cell. 



